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BEAMR - AN INTERACTIVE GRAPHIC COMPUTER PROGRAM FOR DESIGN OF 
CHARGED-PARTICLE BEAM TRANSPORT SYSTEMS 

by Regis F. Leonard and Charles C. Giamati 
Lewis Research Center 
SUMMARY 

A computer program is presented which calculates the characteristics of a 
charged-particle beam as it is transported through a sequence of focusing and bend- 
ing magnets. The maximum dimensions of the beam envelope normal to the trans- 
port system axis are continuously plotted on an oscilloscope as a function of distance 
along the axis . Provision is made to iterate the calculation by changing the types of 
magnets , their positions , and their field strengths . The program is especially 
useful for transport system design studies because of the ease and rapidity of alter- 
ing parameters from panel switches . A typical calculation for a system with eight 
elements is completed in less than 10 seconds. A version containing more-detailed 
printed output but no oscilloscope display is also presented . 

INTRODUCTION 

Designing a beam-handling system for the experimental facility associated with 
the modified NASA Lewis 60-inch cyclotron required that we calculate the effect of 
each proposed beam-handling element upon the beam properties . In fact , calcu- 
lating the properties of the beam at each point in the system from the source (the 
cyclotron) to the target was necessary in order to determine not only the character- 
istics and physical locations of the beam-handling equipment , but also the field- 
settings which would provide the highest quality beam for a particular experiment 
after the hardware had been placed . 

Four characteristics of a charged-particle beam are most significant to a pro- 
spective user . The first is the spatial extent of the beam (x ,y) , both as it traverses 
the beam-handling system and as it arrives at the target. Closely related to the 
spatial extent , and for many experiments of equal importance , is the divergence of 
the beam (x* ,y') when it reaches the target. In addition, it is often necessary to 



know the spread in momentum Ap of the particles contained in the beam . For experi- 
ments which involve timing , it is also necessary to have some knowledge of the 
spread in time (anisochronicity A\) of the incident beam . Thus , a complete des- 
cription of the beam , at least for purposes of any experiment presently anticipated , 
requires six dimensions. 

MATHEMATICAL METHOD 
Single Charged Particle 

The basic assumption underlying all calculations performed by the present 
program is that the transformation of the coordinates of a particle at any point in the 
system to a later point is linear ( i.e. , the calculations are carried out to first order) . 
The transformation may be written 



where X.^ is a six-dimensional vector whose components are the coordinates of a 
particle at time (or position) 1 , is a vector comprised of the set of these same 
coordinates at time (or position) 0 , and is the transformation matrix . The scalar 
coordinates selected for description of a single particle traversing the beam- 
handling system are 

x distance of particle from system centerline in the horizontal plane 
x' angle between system centerline and projection of particle momentum onto 
horizontal plane; x' = dx/dz, where z is coordinate along system axis 
Ap deviation of particle momentum from that of particle which would traverse 
system centerline 

A\ deviation in distance traveled from that of particle traversing system center- 
line 

y distance of particle from system centerline in vertical plane 
y' angle between system centerline and projection of particle momentum onto 
vertical plane; y' is analogous to x' 

The coordinates x, x’ , y , y' are shown graphically in figure 1. 

The problem then becomes one of determining how each type of beam-handling 
element encountered will transform the particle coordinates (i.e. , what is the nature 


2 



of the transformation matrix & for each piece of hardware) . These matrices are 
well known (refs . 1 to 3) and are easily constructed . For the convenience of the 
reader, they are summarized here as follows: 

(1) For a field-free region (drift space) of length L, the & matrix is given by 

1 L 0 0 

0 10 0 
g- _ 0 0 10 

0 0 0 1 

0 0 0 0 

\ 0 0 0 0 

Hence, Xj = is equivalent to the scalar equations 

X 1 = x 0 + Lx 0 

APi=Ap 0 

A\ = A\ 

1 0 

y l = y o + Ly o 

y l =y 0 

(2) For a uniform-field magnet, the & matrix is given by 
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R sin <p 


R(1 - cos cp) 


cos( ip — fi) s j n ^ + tgjj _ cos y) 

cos p 

0 1 
R{1 - cos ip) 


R (<p - sin <p) 
0 
0 


0 0 

1 0 

0 1 - <p tan a 

0 -tan a - tan p + cp tan a tan p 


where R, ( P, a, and /3 are magnet dimensions, defined by figure 2. 

(3) For a quadrupole lens - horizontally defocusing, vertically focusing (ref. 2), 
the matrix is given by 



jL sinh(kL) 0 0 

cosh(kL) 0 0 

0 10 


0 

0 

0 


0 0 
0 0 


0 

0 

0 

0 

cos (kL) 
-k sin(kL) 



where L is the effective length of the lens , and k is the focusing strength related 
to the magnetic field gradient G by 


k = eG/p 


where p/e is the magnetic rigidity . If p/e is expressed in kilogauss-inches , 
and G in kilogauss/inches , then k has the dimensions (in.) 

(4) For a quadrupole lens - horizontally focusing, vertically defocusing (ref. 2), 
the' matrix is given by 
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(5) For a rotation about the beam axis (ref. 3) , the matrix is given by 



where 0 is the angle of rotation - positive for a clockwise (looking downstream) 
rotation, negative for counterclockwise (looking downstream) rotation. The rota- 
tions are mathematical operations in the program which rotate the beam in order to 
use a single bending magnet matrix for all physical orientations of the magnets . 

Beam of Particles 

The formulation just given is sufficient for the calculation of the trajectory of a 
single charged particle through a system of magnets. However, the physical entity 
to be treated is a collection of particles with a spread in initial coordinates . In one 
dimension, for example, the beam would occupy a certain region of the x,x* space. 
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This is most conveniently represented by a phase-space ellipse, as shown in fig- 
ure 3. At a waist, or spatial minimum, the phase-space ellipse becomes upright. 

As the beam traverses a field-free region , for example , the ellipse becomes skewed , 
as shown by the third ellipse in figure 3. 

The formulation of the transformation of a beam envelope or ellipse is again 
most easily accomplished by means of a matrix formalism (ref. 3) . If the beam en- 
velope at a time (or position) 0 is given by 

axx + bxx' + cx'x' = 1 (1) 

then the envelope at time (or position) 1 will be given by a similar ellipse but with 
coefficients a' ,b' ,and c' . In matrix form, equation (1) may be written 



or 

x T y _1 x = i 

It is easily shown that a matrix which transforms the coordinates x,x' to new 
values will transform the matrix Sf as 

5^(2) = ^(1)/ T 

Furthermore , the physical dimensions of the beam are related to the diagonal ele- 
ments of the y matrix: 


x = (S 1t ) 
max 11 


1/2 


*hiax (S 22 5 


1/2 


The generalization of the method for a six-dimensional problem is straightforward. 
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The final transport to the target is accomplished by making the target element 
equivalent to a rotation of zero degrees preceded by the actual drift length. 

The program flow is diagrammed in figure 4. Initially, the input parameters 
are typed in for the following quantities (see the section Input Values for details): 

(1) The initial beam vector, XMAX 

(2) The number of elements, JMAX 

(3) The parameters for each element 

The parameters XMAX and JMAX are read in by the main program . Each set of 

parameters for an element is read in by the subroutine REED; the transport matrix 

for that element is written , if desired , by the subroutine RITE . The matrices are 

calculated by the subroutines MAG, QUAD, and ROT. 

For the transport calculation , the index J is incremented by 1 until JMAX is 

reached . For a given J , the subroutine DRIFT calculates the values of x , y 

max max 

in the drift space . The total transport matrix representing the entire system up 

to the given point , is updated after the return from DRIFT . The transport matrix 

is again updated by MATMLT , which multiplies & by the proper element matrix; and 

these values of the ^matrix (at the end of an element) are saved. 

At the end of the calculation , the DSPLA subroutine displays the results and 

pauses (Pause 2) until a further command is entered from the panel (fig. 5) or until 

the light-pen switch is pressed. The oscilloscope display (fig. 6) includes points 

representing x and y as a function of distance z from the entrance of the 
m&x max 

transport system (see the section Display) . If the light-pen is held over a selected 
point and the pen switch depressed (initiating a call to the subroutine LANDT) , the 
& matrix and the beam vector XMAX at that point will be printed on the Teletype- 
writer . Following Pause 2 , the calculation will be repeated or halted according to 

the setting of the repeat toggle switch . In either case , the values of x and y 

max - max 

will be printed on the Teletypewriter if the print toggle switch is set "up . " The 
values printed are for a range of J as selected on the rotary switches . The sub- 
routine OUTPT accomplishes this printing. Pressing the continue button during 
Pause 2 resumes the processing according to the options set on the panel . 

If the calculation is to be repeated , more input is required . If only the magnetic 
field strengths for two quadrupoles are to be changed , this can be done from the 
panel before the continue button is pressed. For changing other parameters, more 
input through the Teletypewriter is required. A (possibly new) value of JMAX is 
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input , along with the number of the element whose parameters , and possibly whose 
type , are to be changed . After the input for an element , the value of the variable 
AGAIN is set to TRUE if another element is to be changed , and is set to FALSE other- 
wise. In the program, the subroutine REPET controls the input and repeats the 
calculation . 

When the calculation is halted , either the present case can be dumped or a new 
calculation can be initiated . Pushbutton D selects the dump . Pushbutton 7 initiates 
the new calculation , which then requires a new input of the beam parameters . 

If a program has been dumped , the calculations can be resumed at a later time . 
The subroutine RFRESH provides an option for printing out the parameter values 

I 

currently stored in the program when it is being reloaded from a prior dump . The 
value of the variable PRINT is entered as TRUE on the Teletypewriter if the printout 
is wanted, FALSE otherwise. 


Display 

The oscilloscope display (fig. 6) consists of 

(1) Values of x and y as function of distance z, where x is plotted as 

max max max 

a positive value and Y max as a negative value from the centerline (The values are 
displayed at 10 equally spaced points in each drift space.) 

(2) Two horizontal lines, corresponding to the physical beam tube wall, a 
physical distance of 3.81 centimeters (1.5 in.) from the beam axis 

(3) Labels for each element (M,Q,R,T) along with an element index number 

(4) A magnified display of the target at the extreme right of the screen (Two 

fiducial lines each represent a distance of 1 mm from the beam centerline . The 

x and y of the beam at the target are displayed at the magnification corres- 
max max 

ponding to the fiducial lines . ) 

Panel 

The panel (fig. 5) consists of 

(1) Lamps - Lamp 1 being on signals Pause 1 at the beginning of the program. 
The continue button is used to proceed . Lamp 2 being on signals Pause 2 at the 
completion of a new calculation and display . The continue button is used to pro- 
ceed. Lamp 7 being on signals Pause 7 for the selection of the "Dump or go to 1" 
option. Either the dump or go-to-1 button is used to proceed. 
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APPENDIX A 

FORTRAN LISTING FOR PDP-1 5 VERSION 


. . . .MAIN FRCGFAH FOP EEAKR, EEAH TRANSPORT 

LOGIC Al BUTTNS , HALT ,TCGGLS, EITEST 
COHMON/ADZI/A (6 ,|6 , 1 5) ,DIST(15) ,ZL(16) ,IEL 

COHMON/ELEMT/RO (15) , P H I ( 1 5) , AL F ( 1 5) , BET (15) ,EL (15) ,FG ( 15) , 
1ALPH (15) ,IQNR (15) 

COHHON/TAX/T (6/, 6) ,XHAX (6) 

COHHON/SAV/TS/ (6,6,15) ,ZNCT (16) 

CONHON/DSPXY/XDSP (150) , YD SP (150,2) ,JBAX 


/.READY, INITIALIZE ARRAYS 

101 CALL LAHPON(I) 

IBUT=0 / 

110 IF (BOTTNS (IBOT) j GO TO 110 

IF (IBUT.LT. 128) GC TO 110 


CALL ZERO (A) 

CALL ZERO (XDSP) 

CALL ZERO (YDSP) 

CALI LAflPGF ( 1) 

INPOT BEAR ,ELENE!iT PAFAHETEFS 

CALL INLIST (XMAX,5H XRAX) 

CALL INPUT1 (JHAX,4HJHAX,1) 

DO 115 J= 1 , JMAX 
CALL REED (J) 

CALL TOGGLS (ITOG) 

IF(. NOT. BITEST (ITOG, 1C) ) GC TC 115 
CALL RITE (J) 


115 CONTINUE 

C INITIAL VALDES FOR FIRST CALCULATION 

J=1 

200 DO 205 M 1= 1 , 6 

DO 205 M2= 1,6 
T (M 1 , M2) =0 . 

205 IF (H1.EO.B2) T(H1,M2) = 1. 

BEGIN LOOP FCP TRANSPORT 

251 DEL= ZL (J) /10. 

IZ=J*10 * 10 
ZTOT = ZNOT(J) 

IF(J.EQ.I) ZTOT = 0. 

252 CALL DRIFT (DEL, ZTOT, IZ) 
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DO 255 L=1,6 

T(1,L) = T(1,L) + T (2 , L) *ZL ( J) 

T (5,L) = T(5,L) ♦ T (6 , I) *ZL { J) 

ZTOT = ZTCT + ZL(J) 

CALL BATHLT(J) 

ZTOT = ZTCT + DIST (J) 

JP 1 = J+1 
ZNCT(JPI) = ZTOT 
DO 275 H 1=1,6 
DO 275 B2= 1 , 6 
275 TSV(B1,B2,J) = T(B1,B2) 

J = JP 1 

29 S IF (J . LE. JBAX) GO TO 251 

C END TRANSPORT LOOP 

301 CALL DSPLA 

CALL TOGGLS (ITOG) 

IF(. NOT. BITEST (ITOG, 14) ) GC TC 402 
40C CALL OUTPT 

402 CALL TOGGLS (ITOG) 

HALT= .TRUE. 

IF (BITEST (ITOG, 12) ) HALT=. FALSE. 

IF (HALT) GC TO 599 
500 CALL REPET (JHIN) 

, J=JBIN 

CALL TOGGLS (ITOG) 

IF(. NOT. BITEST (ITOG, 15)) GO TO 504 
CALL CLKOFF 

504 IF (J.EQ. 1) GO TO 200 

JM1 = J-1 

DO 505 fl 1= 1 , 6 
DO 505 B 2= 1 , 6 

505 T (B 1 , B2) = TSV (B1,B2, JB1) 

GO TO 251 

59 S IE=0 

CALL LABPON (7) 

60 C IF (BOTTNS (IP) ) GO TO 600 

IF(IB.NE. 16, AND. IB. NE. 32) GO TO 600 
CALL LABPCF (7) 

IF ( IB • EQ . 16) GO TO 601 
GO TO 101 

601 CALL VPSTOP 

CALL QDONF 
CALL VPSTFT 
CALL PFRESH 
GO TO 301 
END 
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COMMON/ADZI/A (6,6, 15) ,EIST(15) ,ZL(16) ,IEL 

COMMON/ELEMT/RO (15) ,PHI(15) ,ALF(15) , EET (15) ,EL(15) ,FG(15) , 
1AIFR (15) ,IQMR (15) 

C 

JW = J 

IEL = IQMR(JW) 

WRITE (6,10) JW,IEL,ZL (J) , DIST (J) 

10 FORMAT (1H0. 17H ELEMENT NUMBER =,I1,7H TYPE= ,I1,12H 

PRE-DRIFT= , 1F10. 5, 1 1H LENGTH =,E10. 5/) 

DO 1 i5 K=1,6 

WRITE (6,20) (A (K, N, J) , N= 1 , 6) 

115 CONTINUE 

RETURN 

20 FORMAT ( 1H , F 1 1 . 6 , 5F 1 2 . 6) 

END 


SUBROUTINE MAG (J, BC, PHI, ALP, EET) 

.... COMPOTE MAGNET MATRIX FROM INPUT PARAMETERS 
COMMON/ADZI/A (6,6,15) ,DIST(15) ,ZL(16) ,IEL 

C SET ELEMENTS USEE IN QUAD OR ROT TO ZERO 

A ( 1, 5, J) = 0. 

A (2 , 6 , J) = 0. 

A(5,1,J) = 0. 

A (6 , 2 , J) = 0. 

C . . . . COMPUTE SINE AND CCSINE VALUES 

RPI=3 . 141593/180. 

PR=PHI*RPI 
AR=ALP*RPI 
BP=BET*HPI 
CP=C0S (PR) 

CA=COS ( AR) 

CB=COS (BR) 

SP=SIN (PR) 

SA=SIN (AR) 

SB= SIN (BR) 

TA=SA/CA 

TB=SB/CB 

A(1,1,J) * CF+SP*T A 
A ( 1 , 2, J) = RO*SP 
A (1 , 3, J) = RO* ( 1 . -CP) 

A (2, 1 , J) = -SIN (PR-AR-BR) /(RC*CA*CE) 

A (2 ,2 ,J) = CP ♦ S P*T E 
A(2,3,J) = SP + ( 1 . -CP) *TB 
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A(3,3,J) 
A(4, 1,J) 

a (4,3, J) 
a (4, 4 , j) 
a <5,5 ,j) 
a (5,6,J) 
a (6,5,j) 
a (6,6,J) 
DIST (J) 
RETURN 
END 


1 . 

SP + (1.-CF)*TA 
RO* ( 1 , -CP) 

RO* (PR-SP) 

1 . 

1.-FR*TA 

PR*RO 

(-TA-TE+PR*TA*TE) /RO 
1.-PR*TB 
a (5,6, J) 


SUBROUTINE QO AE (J , EL , F6) 

. . . • COHPUT ooar MATRIX from input parameters 
COM MON /ADZ I /A (6,6,15) , BI ST ( 1 5) ,ZL ( 1 6) ,IEL 

C SET ELEMENTS USED IN MAG OR ROT TO ZERO 

A ( 1 , 3 , J) = 0. 

A (1 ,5, J) = 0. 

A (2, 3, J) = C. 

A (2 ,6 , J) = 0. 

A (4, 1 , J) = 0. 

A (4 ,2, J) = 0. 

A (4, 3, J) = 0. 

A (5 , 1 , J) =0. 

A (6, 2,0) = 0. 

. CGHPfJTE MATRIX ELEMENTS 

IF(FG) 100,200,300 
100 FG = - FG 

EPL = EXP (FG*EL) 

EMI = EXP (-FG* EL) 

SF = SIN (FG*EL) 

A ( 1 , 1 , J) = .5* (EPL + EMI) 

A (1 , 2, J) = (EPL-EMI) / (2.*FG) 

A (2, 1 , J) = FG*FG*A ( 1, 2, J) 

A (5,5, J) = COS (FG*EL) 

A (5,6, J) = SF/FG 
A (6,5,0) = -FG*SF 
FG = -FG 
GO TO 400 

30 C SF = SIN (FG*EL) 

EPL = EXP (FG*EL) 

EMI = EXP (-FG* EL) 
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A (1, 1,J) = COS (FG*El) 
1(1,2,J) = SF/FG 
A (2,1,0) = -FG*SF 
A(5,5,J) = .5+ (EPL+E N I) 

A (5*6 , J) = (EPL-EMI) / (2.*FG) 
A (6* 5* J) = FG*FG*A (5,6, J) 

GO TO 400 

20C A (1, 1, J) = 1. 

A (1,2,0) = El 
A (2, 1 , J) = 0. 

A(5,5,J) = 1. 

A (5, 6 , J) = EL 
A (6 , 5 , J) = 0. 

400 A (2, 2, J) = A ( 1 , 1 , J) 

A(3,3,J) = 1. 

A (4, 4, J) = 1. 

A (6 ,6 , J) = A (5,5, J) 

DIST(J) = EL 

RETURN 

END 


SUBROUTINE HOT(J,AA) 

. . . .COMPUTE ROTATION MATRIX FROM INPUT FARAMATERS 
COMMON/ ADZI/A (6,6, 15) , EIST ( 1 5) , ZL ( 16) ,IEL 

C SET ELEMENTS USED IN MAG OR QUAE TC ZERO 

A (1 ,2, J) = 0. 

A ( 1 , 3, J) = 0. 

A (2 , 1 , J) = 0. 

A (2, 3, J) = 0. 

A (4 , 1 , J) = 0. 

A (4, 2, J) =0. 

A (4,3, J) = 0. 

A ( 5 , 6 , J) =0. 

A (6 ,5 , J) = 0. 


C COMPUTE MATRIXEIEMENTS 

AP = AA*3. 141593/180. 

CS = CCS (AP) 

SN = SIN (AR) 


A ( 1 , 1 , J) = CS 
A (1 ,5, J) = SN 
A (2, 2, J) = CS 
A (2 ,6 , J) = SN 
A (5, 1 , J) = - SN 
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A (5,5, J) = CS 
A (6,2,3) « -SB 
A(6,6,J) = CS 
A (4, 4, J) = 1. 

A (3 ,3 , J) = 1. 
DIST(J) = 0. 
RETURN 
END 


SUBROUTINE DRIFT (DEL,7TCT,IZ) 

. . . .TRANSPORT BEAN THRU DRIFT TUEE, CONFUTE ECSITIONS 
DIMENSION SIG1 (6,6) ,SIG2 (6,6) ,XMEX (6) 

COHMON/ADZI/A (6 , 6 , 1 5) , DI ST ( 1 5) , ZL ( 16) ,IEL 
COHNON/DSFXY/XDSP ( 150) , YDSP ( 150, 2) , JM AX 
COHHON/TAX/T (6,6) , X MAX (6) 

COMMON/S AV/TSV (6,6,15) ,ZNCT(16) 

EQUIVALENCE (A (1,1,15) ,SIG1 (1,1) ) 

EQUIVALENCE (TSV (1, 1, 15) ,SIG2 (1, 1)) 

ZINT=ZTOT 
DO 215 M 1= 1 , 6 
DO 215 M 2= 1 , 6 
SIG 1 (N 1 , M2) =0 . 

DO 205 M3= 1 , 6 

SIG 1 (Hl,H2)=SIGl (HI, M2) +T (M 1 , M3) * (XMAX (M3) **2) *T (M2, M3) 
205 CONTINUE 

SIG2 (H1,M2)=SIG1 (Ml, M2) 

215 CONTINUE 

DO 255 K 1= 1 , 10 
DO 225 M 1= 1 , 6 

SIG2 (1,H1)=SIG2 (1,M1) 4DEL*SIG 1 (2, M 1) 

SIG2 ( 5, M 1) =SIG2 (5, Ml) 4DEL*SIG1 (6, Ml) 

225 CONTINUE 

DO 235 M 1= 1 , 6 

SIG2 (HI, 1)=SIG2 (Ml, 1) +DEL+SIG 1 (M 1,2) 

SIG2 (M 1 , 5) =SIG2 (Ml ,5) 4DEL*SIG1 (Ml ,6) 

235 CONTINUE 

SIG2 (1, 1) =SIG2 (1, 1) + (DEL**2) * SIG 1 (2,2) 

SIG2 (1,5)=SIG2 (1,5) 4 <EEL**2) *SIG 1 (2,6) 

SIG2 (5,1) =SIG2 (5,1) + (DEL** 2) *SIG 1 (6,2) 

SIG2 (5,5) =SIG2 (5, 5) 4 (DEL**2) * SIG 1 (6,6) 

IF (IZ • EQ. 999) GO TO 260 

ZINT-ZINT4DEL 

IZ=IZ4l 


18 



YDSP(IZ,1) = SQRT (SIG2 (1,1) > 

XDSP (IZ) = ZINT 

YDSP ( IZ , 2) = SQRT (SIG2 (5,5) ) 

DO 245 M1=1,6 
DO 245 H2=1,6 
SIG 1 (M1,M2)=SIG2(M1,M2) 

245 CONTINUE 

255 CONTINUE 

RETURN 

26C DO 265 ID=1,6 

265 XMEX(ID) = SQRT (SIG2 (ID, ID) ) 

WRITE (6,10) (XHEX(IW) ,IW=1,6) 

10 FORMAT (1H0,6HXMAX= ,6(3X,F7 .3)) 

RETURN 
END 


SUBROUTINE MATMLT (JJ) 

C 

C . . . .ROUTINE FOR MATRIX MULTI PLICATION 
DIMENSION C (6 , 6) 

COMMON/ADZI/A (6 , 6 , 1 5) , DI ST (1 5) , Z1 ( 16) ,IEL 
COMMON/TAX/T (6,6) ,XMAX (6) 

EQUIVALENCE (A ( 1 , 1 , 1 5) ,C (1 , 1 ) ) 

DC 105 M1= 1 ,6 
DO 105 M2= 1 , 6 
C (M 1 , M2) =0 . 

DO 105 M 3= 1 , 6 

C (M1,M2) = A (HI, M3, JJ) *T (M3, M2) ♦C(H1,H2) 

105 CONTINUE 

DO 115 M 1= 1 ,6 
DO 115 M2= 1 , 6 
T (M 1 , M2) =C (Ml, M2) 

115 CONTINUE 

RETURN 
END 


SUBROUTINE DSPLA 
C 

C THIS ROUTINE USES SCOPE DISPLAY 

INTEGER XYP (2,320) ,YPX,YPY,V (5,2) ,P (2, 2) , TAGS (4, 30) 
LOGICAL BUTINS, PENSH 
DIMENSION RJ ( 15) ,RL(4) 
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n o 


COMMON/TAG/RT AGS (2,30) ,XYP 
COMBON/DSPXY/XDSP (150) ,YDSP{150,2) ,JMAX 

COBBON/ELEBT/BO (15) ,PHI (15) , ALP (15) ,EET (15) ,EL(15) ,FG(15) , 
1ALPH (15) ,IQMR (15) 

EQUIVALENCE (TAGS (1,1) , STAGS (1,1)) 

DATA RL (1) ,PL (2) , RL (3) , PL (4) / IBM , 1HQ, 1 ER , 1HT/ 

DATA RJ ( 1) ,F«J (2) ,RJ (3) ,RJ(4) ,BJ(5) ,PJ<6) ,RJ(7) ,P J (8) ,RJ(9) , 
1RJ (10) ,RJ(11) , RJ (12) , R J (13) , P J (14) , R J (15) /1H1,1H2,1H3,1H4, 
21H5, 1R6, 1H7, 1H8, 1H9 , 2H 10 , 2H 1 1 , 2B 1 2, 2H 13, 2H14,2H15/ 

Z SCALE, 1 FOINT/INCH 

VERTICAL SCALE, 1 PCINT=0.005 INCH (1.5"=300) 

ICC NP= JBAX* 1 0 

DO 125 K= 1 , NP 

KY=2*K 

KX=KY-1 

XYP (2, KX) =512.+200.*YDSP (K,1) 

XYP (2,KY)=512.-200.*YDSP (K,2) 

XPZ=XDSP (K) 

XYP (1 ,KY) = XPZ 
125 XYP ( 1 , KX) = XPZ 

NP2=2*NP 
NX=NP2 ♦ 1 
NY=NX ♦ 1 

XYP (2, NX) = 512. + 1000.* YD SP (NP,1) 

XYP (2 , NY ) - 512. - 100C.*YDSP (NP,2) 

YPX = XYP (2, NX) 

YPY = XYP ( 2 , N Y) 

XYP (1, NX) = 1000 
XYP ( 1 , N Y) = 1022 

C GENERATE 1 BB BARKERS 

20C NB 1 = NY+1 

NB2 = NB 1 ♦ 9 
MB = 40 
NZ * 1010 
DO 225 L=NB1,NB2 
XYP (2 , L) = 512 + MM 
XYP ( 1 , L) = NZ 
MM = - MM 

225 NZ = NZ ♦ 1 

NP2 = NP2 12 

CALL VPPNTS(XYP,NP2,IHIT,1) 

C GENERATE ECUNDARIES AT ♦ 1 . 5" , - 1 . 5" 

300 MINUS = -1 

DO 325 LL= 1,2 
P(1,1) = 0 
P (1 , 2) = 950 
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HP = 512 + MINUS*300 
P (2, 1) = HP 
P (2, 2) = HP 
CALL VECTOR (P r V r LL) 

325 MINOS = - MINOS 

39 C CALL VPVECT (V,2,IH,2) 

4CC MO = 1 

DO 455 J= 1 # JMAX 
JE = 2*J 

JO = JE- 1 

RTAGS(1,JO) = RJ(J) 

K = (IOHP (J) + 1) / 2 
PTAGS(1,JE) = RL (K) 

4 1 C J10 = 10*J 

TAGS (3, JO) = XDSP (J10) 

TAGS (3 r JE) = XDSP (J 10) 

IF(IQMR (J) .NE. 3) GO TC 420 
MQLO = 0 

IF(MQ.EQ.-I) MQLO = 3C 

TAGS (4, JO) = 512 ♦ MQ*(305 + HOLO) 

TAGS (4, JE) = 512 ♦ MQ* (305 + 30 ♦ MQLO) 
MQ = - MO 
GO TO 455 

420 TAGS (4, JO) = 512 + 305 

TAGS (4, JE) = 512 + 305 + 30 
455 CONTINOE 

NT = JMAX * 2 

49 C CALL VPTAGS (TAGS , NT, IT, 3) 

500 CALL VPSTPT 

CALL CLKON ' 

IF (IPX. GT. 552. OR .YPY. IT. 472) GO TO 526 
DO 525 LA=1,7 

525 CALL LAMPON (LA) 

526 CALL LAMPCN (2) 

IBOT=0 

53C IF (PE NSW (.TROE.) ) GC TO 600 

IF(BUTTNS (IBUT) ) GO TC 530 
IF (IBOT.LT. 128) GO TO 530 
DO 535 LA=1,7 
535 CALL LAMPCF (LA) 

BET URN 

600 NH= (IHIT+ 1 ) /2 

CALL LANDT(NH) 

GO TO 530 
END 
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SUBROUTINE COTPT 


C 

C . . . . BRITE COT DISPLAY POINTS 

COMMON /DSPXY/XDSP (150) ,YDSF(150,2) , JMAX 
BRITE (6,10) 

10 FORMAT (1HG,5X,45HX AND Y BEAMCCCFDSAS FUNCTION OF DISTANCE Z/) 

J3 = IFOTOE (3) 

J 4 = IFOTOF (4) 

JBGN = 10*03 - 9 
JEND = 10* J4 

IF (JEND.GT. 140) JENB=140 
DO 115 J1=JBGN,JEND, 1C 
J2= J1 *9 

BRITE (6, 20) (YDSP(JX,1) ,JX=J1,J2) 

20 FORMAT (2H X,10F7.3) 

WRITE (6,30) (XDSP(JZ) ,JZ=J1,J2) 

30 FORMAT ( 1H ,10F7.2) 

BR ITE (6 , 40) (YDSP (JY ,2) ,JY=J1,J2) 

40 FORMAT (2H Y,10F7.3///) 

115 CONTINUE 

RETURN 
END 


SUBROUTINE REPET (JMIN) 

C 

C . . . .CHANGES SPECIFIED ELEMENTS AND PEPEATS TRANSPORT 
LOGICAL AGAIN, TOGGLS, EITEST 

COMMON/ELEMT/PO (15) ,PHI (15) , ALP (15) ,PET (15) ,EL (15) , FG ( 15). , 
1ALPH (15) ,ICMR (15) 

COMMON/DSP XY/XDSP (150) , YDS F (150,2) ,JMAX 
C 

C IF TOGGL 11 IS SET, GET NEB QUAD CURRENT FROM THUMB 

CALL TOGGLS (ITOG) 

IF(. NOT. BITEST (ITOG, 1 1) ) GC TG 301 
J1=IROTOR (1) 

J2=IFOTOR (2) 

JMIN = J1 

IF(J2.LT.J1) JMIN = J2 
IFG 1 = ITHUMB (1) 

IFG2=ITHUMB (2) 

FG1=FIOAT (IFG 1) *0.0000 1 
FG2=FLOAT (IFG2) *0 . 000C 1 
IF (BITEST (ITOG ,13)) FG1=-FG1 
IF (BITEST (ITOG, 17) ) FG2=-FG2 
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10 


FG (J 1 ) = FG1 
FG (J2) = FG2 
El 1 = EL { J 1 ) 

EL2=EL (J 2) 

WRITE (6,10) J1,FG1, J2,FG2 
FORMAT (1HC,6HJ, FG,5X,2(I2,F10.5)) 
CALL QUAD (Jl, ELI, FG1) 

CALL QUAD (J2,EL2,FG2) 

IF (.NOT. BITEST (ITOG, 10)) GO TO 50C 
CALL BITE (J 1) 

CALL BITE (J2) 

GO TO 500 
C 


C INPOT ELEMENT PARAMETERS THBO TELETYPE 

301 CALI INPUT1 (JMAX,4HJMAX, 1) 

JMIN = JMAX 

302 CALL INPOT1 (JJ,5HNELEM, 1) 

IF (JJ.LT. JMIN) JMIN=JJ 
CALL BE ED (JO) 

CALL TOGGLS (ITOG) 


IF(. NOT. BITEST (ITOG, 10)) GO TO 303 
CALL BITE (JJ) 

303 CALL INPOT1 (AGAIN, 5BAGAIN, 4) 

IF (AGAIN) GO TO 302 
500 FETUBN 

END 


SOBBOOTINE RFRESH 
C 

C WRITE PARAMETERS AFTER RESTART FROM PREVIOUS DUMP 

LOGICAI FFINP 

COMMON/ELEMT/RO (15) ,PHI (15) , ALE (15) , BET (1 5) , EL ( 15) , FG ( 15) , 
1 ALPH (15) , IOMR (15) 

COMMON/ADZI/A (6,6,15) ,DIST(15) ,ZL(16) ,IEL 
COMMON /DSPXY/XDSP (150) ,YDSP (150,2) ,JMAX 
C 

WRITE (6,20) 

20 FORMAT (42H0 DO YOU WISH TO PRINT CURRENT PARAMETERS?) 

CALL INPUT1 (PRINP,5HPRINT,4) 

IF (.NOT. PRINP) GO TO 500 

DO 295 L= 1 , JMAX 

LW=L 

WRITE (6,40) LW,IQMR (L) ,ZL(L) ,EIST(L) 

40 FORMAT (1H0, 17H ELEMENT NUMBER =,I2,6H TYPE ,I1,12H 

PRE-DBIFT= , 
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1F10.2,11H LENGTH- , F 10 . 2) 

IW=IQMR(L) 

GO TO (231,231,233,233,235,235,5C0),IW 
231 SHITE (6, 10) RO (L) ,PHI (I) ,AIF (L) ,BET (L) 

10 FORMAT (1H0,6H RHO = , F 10 . 2 , 6H PHI = , F 1 0 . 2, 8H ALPHA= , F 1 0. 2 , 7H BE 

1TA =,F10.2/t 
GO TO 295 

233 WRITE (6, 30) FG (L) 

30 FORMAT (1H0, 11H FLD GRAB =,F10.5/) 

GO TO 295 

235 WRITE (6,50) ALPH(L) 

50 FORMAT(1H0,17H ROTATION ANGLE =,F10.2/) 

295 CONTINUE 

500 RETURN 

END 


SUBROUTINE LANDT(NH) 

C 

C CALCniATE I FROM XMAX; CALCULATE T MATRIX 

COMMON /ADZI/ A (6,6,15) ,DIST(15) ,ZL(16) ,IEL 
COMMON/TAX/T (6,6) ,XKAX (6) 

COMMON/ELEMT/RO (15) ,PHI(15) , ALP (15) , BET (15) , EL ( 1 5) , FG ( 1 5) , 
1ALPH (15) , IQMR (15) 

COMMON/S A V/TSV (6,6,15) ,ZNCTt16) 


C JH IS INDEX FOR HIT 

100 JH= (NH-1) /10 + 1 

JHM1 = JH - 1 
DEL=ZL (JH) /10. 

IZL= NH- JHM 1* 1 0 
ZLH = CEL*FLCAT (IZL) 

IF ( JH . NE . 1) GO TO 250 

C .SET UP SPECIAL CASE FOP JH=1 

JHM 1= 1 5 

210 DO 215 M 1= 1 , 6 

DO 215 M2=1,6 
TSV (Hi ,M2,15) = 0. 

215 IF (Ml. EQ .M2) TSV (W 1 , M2 , 15) = 1. 

C STORE TRANSPORT MATRIX IN MI, ft, 15) 

250 DO 255 L= 1 , 6 

A (1»L,15) = TSV ( 1 , L , JEM 1) ♦ TSV (2, L, JHM 1) *ZLH 
A(2,L,15) = TSV (2 , L , JHH1) 

A(3,L,15) = TSV (3, L, JHM 1) 

A (4, L , 15) = TSV (4,L, JHM11 

A (5, L, 15) = TSV (5 , L , JHM 1) ♦ TSV (6, L , JHM 1) *ZLH 
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255 A(6,L,15) = TS V (6 , L, JHM1) 

C SET OP PARAMETERS EOF CALL TO RITE 

DIST (15) = ZLH 
IQMR(15) = 9 
ZL ( 15) = ZLH 

IP(JH.NE.I) ZL (15) = ZL (15) + ZNOT (JR) 

CALL RITE (15) 

C SET UP T MATRIX FCR CALL TO DRIFT 

DO 305 M1=1,6 
DO 305 M2= 1 ,6 

305 T(M1,M2) = TSV (M1,W2, JHM1) 

DEL=ZLH 

CALL DRIFT (DEL, 0. ,999) 

RETURN 

END 
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APPENDIX B 


IBM 7094 VERSION OF BEAMR 

The IBM 7094 version of the BEAMR program is similar to that just described 
except that there exists no interactive mode and no facility for display of results . 
Rather, the program supplies printed description of the beam and cumulative values 
for the transfer matrix at points specified. In addition, the program will auto- 
matically search to determine what values of the quadrupole focusing strength are 
required to obtain preset values for two elements of the total transfer matrix at a 
preset position . 

The input for the program is as follows: 


CARD 

NUMBER 

FORMAT 

INFORMATION 

1 

13A6 

Title 

2 

6F10.0 

^ax’^nax’^^’^’^max’ ^max 

3 

212 

JMIX , JMOX 

4 

1F10.0 

DELK 

5 

5F10.0.I1 

P1,P2,P3,P4,P5,XJ 

4+JMIX 

- 

- 

4+JMIX+l 

413 , 4E 15 . 8 , 213 AI ,AJ ,BI ,BJ , SFOC 1 , SFOC 2 , VAL1 , VAL2 , J 1 , J2 

(This card only for JMOX=l.) 


Definitions of Input 

The parameters x , x' , Ap/p, A\, y m „ , y* are characteristics of the 

max max max max 

beam at the entrance to the system, as defined previously , in inches and radians. 
JMIX is the total number of beam-handling elements in the system , including drift 
spaces. JMOX specifies search procedures: 0, no search; 1, search on quadrupole 
strengths . 

The total transfer matrix and the maximum beam dimensions in the drift space 
beyond the last beam-handling element will be calculated and printed out 25 times 


26 



if so desired . These calculations will be done at intervals DELK (in . ) from the exit 
of this element . 

The definitions of PI, P2, P3, P4, and P5 vary depending on the type of beam- 
handling element involved. The type of element is defined by XJ: 

For a magnet , XJ=1: 

PI radius of curvature 

P2 angle of deflection 

P3 entrance edge angle 

P4 exit edge angle 

P5 PL 

For a drift space , XJ=2: 

PI length of drift space 

P2 0 

P3 0 

P4 0 

P5 0 

For a quadrupole magnet, XJ=3,: 

PI effective length of quadrupole 

P2 strength, k - positive for horizontal focusing, negative for vertical focusing 
P3 PL 

P4 0 

P5 0 

For a quadrupole (part of a symmetric triplet) , XJ=4: P1,P2,P3,P4, and P5 are the 
same as for XJ=3. For a rotation, XJ=5: PI is the magnitude and direction of the 
rotation - positive for a clockwise rotation , negative for a counterclockwise rotation . 

The quantity PL defines , for each beam-handling element , whether a des- 
cription of the beam for the drift space following the element will be printed out . If 
PL = 1 , the drift space following the element will be divided into 20 equal intervals 
and a description of the beam (total S matrix and maximum beam dimensions) will 
be printed as calculated at the end of each interval . If PL = 0 , no output is printed . 

The quantities AI,AJ,BI,BJ,SF0C1,SF0C2,VAL1,VAL2,J1, and J2 specifiy 
search procedures . The search routine will attempt to adjust the strengths of quad- 
ruples J1 and J2 so as to produce, at a distance SFOC1 beyond the last element, the 
value VAL1 for the transfer matrix element T(AI,AJ) . Simultaneously, it will attempt 
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to produce at the distance SFOC2 beyond the last element the value VAL2 for the 
transfer matrix element T (BI ,BJ) • 

If one of the quadrupoles specified is the first section of a symmetric triplet 
(XJ=4) , the focusing strength of the last section of the triplet will be held identically 
equal to that of the first section during all search procedures . 

The search routine optimizes the specified parameters by first constructing 
numerically the gradients 


9T(AI, AJ) 
9k(Jl) 


9T(AI, AJ) 
9k(J2) 


9T(BI, BJ) 
9k(Jl) 


9T(BI, BJ) 
9k(J2) 


and then solving the equations 


AT(AI, AJ) = 

9T(AI, AJ) 

Ak(Jl) + 

9T(AI, AJ) 

Ak(J2) 

9k( J1 ) 

9k(J2) 

AT(BI, BJ) = 

» — 

9T(BI, BJ) 

9k(Jl) 

Ak(Jl) + 

9T(BI,BJ)1 

9k(J2) 

Ak(J2) 


for corrections Ak to the field gradients. The quantity AT (AI,AJ) is the difference 
between the current value of the (AI,AJ) matrix element and the desired value VAL1. 
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APPENDIX C 


FORTRAN LISTING OF IBM 7094 VERSION 


$1 EFTC BEAMR DECK 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


BEAMR DOES EEAM TRANSPORT CALCULATIONS IN FIRST ORDER, USING 
A BEAM TRANSFER MATRIX. IT TREATS MAGNETS , DRIFT SPACES , QUACROPOLE 
LENSES, AND ROTATIONS. 

INPUT IS AS FOLLOWS 
CARD ONE CONTAINS A TITLE 

CARD 2 (6F10.0) CONTAINS A DESCRIPTION OF THE BEAM AT THE SYSTEM 
ENTRANCE 

A) MAXIMUM SIZE OF BEAM IN THE X-DIRECTION (INCHES) 

B) MAXIMUM DIVERGENCE OF THE BEAM IN THE X PLANE (RADIANS) 

C) INCIDENT MOMENTUM SPREAD (FRACTIONAL) 

D) MAXIMUM SPREAD IN PATH LENGTH ENTERING THE SYSTEM (INCHES) 

E) MAXIMUM EXTENSION OF THE EEAM IN THE Y DIRECTION (INCHES) 

F) MAXIMUM DIVERGENCE OF THE BEAM IN THE Y PLANE (RADIANS) 

CARD 3 (FORMAT 212) CONTAINS JMIX AND JMOX 

JMIX IS THE TOTAL NUMBER OF ELEMENTS (DRIFT SPACES, QUADRUPOLE 
SINGLETS, MAGNETS, ROTATIONS) IN THE EEAM LINE 
JMOX SPECIFIES SEARCH PROCEDURES 
JMOX=0 IMPLIES NO SEARCH 
JMOX=1 IMPLIES SEARCHING ON QUADRUPOLE CURRENTS 
CARD 4 (1F10.C) CONTAINS THE QUANTITY DELK (INCHES) 

CARDS 5 TO 5* JMIX (5F10. 0,111) EESCPIEE INDIVIDUAL BEAM ELEMENTS, 
ONE CARD PER ELEMENT 

MAGNET- COL 1-10 CONTAIN THE RADIUS OF CURVATURE (GREATER THAN 

ZERO) 

COL 11-20 CONTAIN THE ANGLE OF DEFLECTION (GREATER THAN 

ZERO) 

COL 21-30 CONTAIN TEE ENTFANCE ANGLE (POSITIVE ANGLES 

YIELD VERTICAL FOCUSING) 

COL 31-40 CONTAIN THE EXIT ANGIE ( POSITIVE ANGLES 

YIELD VERTICAL FOCUSING) 

COL 41-50 CONTAIN PL 

COL 51 CONTAINS XJ (SHOUID EQUAL 1 FOR A MAGNET) 

THE AEOVE SIGN CONVENTIONS ARE FOR A MAGNET WHICH 
PENDS TO THE RIGHT. FOR A LEFT-BENDING MAGNET ALL 
SIGNS ARE REVERSED. 

QUADRUPOLE COL 1-10 EFFECTIVE LENGTH OF THE QUAD (INCHES) 

COL 11-20 CONTAIN K , WHERE K IS THE QUADRUPOLE 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


STRENGTH DEFINED EY 

K = SQRT (FIELD GRADIENT (KG/INCH) /BP (KG-INCHES) ) 
POSITIVE K YIELDS X-FOCUS, Y-DEFOCUS 
NEGATIVE K_ YIELDS X-DEFOCUS, Y-FOCUS 
COL 21-30 CONTAIN EL 

COL 51 CONTAINS XU (SRCOLD BE 3 POR A QOADRUPOLE) 

IF COL 51 CONTAINS A 4 IT INDICATES THAT THE ELEMENT 
BEING DESCRIBED IS TO BE CCNSI EERED AS THE FIRST 
SECTION OF A SYMMETRIC TRIPLET, AND IN ANY SEARCH 
PROCEDURE THE LAST ELEMENT OF THE TRIPLET WILL ALWAYS 
HAVE A STRENGTH EQUAL TO THIS ONE 
DRIFT SPACE COL 1-1C LENGTH OF THE DRIFT (INCHES) 

COL 11-21 PL (SHOULD EE ZERO) 

COL 51 XU (SHOULD BE 2 FOR A DRIFT SPACE) 

ROTATION, COL 1-10 CONTAIN THE MAGNITUDE OF THE ROTATION 

IN DEGREES. . A POSITIVE ROTATION INDICATES A 
DOWNWARD EEND, A NEGATIVE ROTATION INDICATES 
AN UPWARD E ENT 
COL 17-20 CONTAIN PL 

COL 51 CONTAINS XJ (SHOULD BE 5 FOR A ROTATION) 

FOR ANY ELEMENT, SETTING PL=1 CAUSES THE TRANSFER MATRIX TO EE 
PRINTED CUT AT INTERVALS OF DELTA INCHES BEYOND THE GIVEN 
ELEMENT, WHERE DELTA=1/20 OF THE FOLLOWING DRIFT SPACE. FOR THE 
LAST ELEMENT THE TRANSFER MATRIX IS PRINTED OUT 25 TIMES IN 
INTERVALS OF DELK. 

IF PL=0 , THE TRANSFER MATRIX IS NOT PRINTED AT THE EXIT CF 
THE GIVEN ELEMENT. 

LAST CARE (413 , 4E 1 5 . 8 , 21 3) 

CONTAINS AI ,AJ,BI,BJ, SFOC 1,SBOC2,VALl ,VAL2,J1,J2. 

THE SEARCH ROUTINE WILL ATTEMPT TO VARY THE STRENGTHS OF THE 
J1-TH AND J2-TH QUADRUFOLES IN CPDER TO MAKE THE (AI,AU) 

AND (PI , EJ) ELEMENTS GF THE TRANSFER MATRIX EQUAL TO THE VALUES 
VAL 1 AND VAL2 RESPECTIVELY AT DISTANCES SFOC 1 AND SFOC2 
RESPECTIVELY BEYOND THE LAST BEAM ELEMENT. 


INTEGER X , F 

DIMENSION X (50) ,RO (50) ,PHI (50) , ALP (50) , BET (50) , ALPP(50) , PHIR (50) , - 
1BETR (50) ,A (6,6,50) ,Z (50) ,T (6,6) , TITLE (13) ,BINT1 (3) ,BINT2(3) , 

2?G (50) , EL (50) , GRAD (2,2) ,TINT1 (6,6,3) , TINT 2 (6, 6, 3) ,PL (50) ,ALPH(50) ,- 
3 ALPHR (50) ,A1 (6,6) ,A2 (6,6) ,XMAX (6) ,AC (5,50) 

COMMON A ,T, XMAX 


100 READ (5,1) 
READ ( 5, 2) 
READ (5,3) 
R EAD (5,2) 


TITLE 

(XMAX (I) ,1=1,6) 
JMIX, JMOX 
DELK 
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Receding page blank not film® 


ALPHR (J ) = ALPH (J) *3.14155265/180. 
A (1 , 1 , J) =CCS (ALPHR (J) ) 

A (1,5, J)=SIN (ALPHR (J) ) 

A (2,2,J) =COS (ALPHR (J) ) 

A (2,6, J) = SI N (ALPHR (J) ) 

A (5, 1,0) =-SIN (ALPHR (J) ) 

A (5,5, J) =COS (ALPHR (J) ) 

A (6,2, J) =-SIN (ALPHR (J) ) 

A (6,6, J)=COS (ALPHR (J) ) 

A (4, 4, J) =1.0 
A (3,3, J) =1 . 0 
WRITE (6, 15) J , ALPH ( J) 


199 WRITE(6,10) 
WRITE (6,10) 
WRITE (6, 10) 
WRITE (6,10) 
WRITE (6, 10) 
WRITE (6,10) 

200 CONTINUE 


(A (1 , N, J) , N= 1 , 6) 
(A (2, N, J) , N= 1 , 6) 
(A (3,N, J) , N= 1 , 6) 
(A (4 , N, J) , N= 1 , 6) 
(A (5,N, J) , N= 1 , 6) 
(A (6 , N, J) , N= 1 , 6) 


203 DO 201 M 1=1 ,6 
DO 201 M2= 1 ,6 
T (Hi, M2) = A (ft 1 , M2, 1) 

201 CONTINUE 

DO 300 J=2,JMIX 
CALL MATMLT(J) 

IF (X (OHAX) . NE.O) GO TO 300 
IF (PL ( J) . EQ . 0. 0) GO TO 300 
WRITE (6,14) J,X(J) 

DEL=(Z (J+1) /20.0) +1.0 

IDEL=DEL 

DEL=IDEL 

IF (J+1.EQ. JMAX) £EL=DELK 
CALL EXIT (DEL) 

300 CONTINUE 


401 IF (X (JM AX) ) 100,100,998 

998 IF (F-2) 999,1000,1000 

999 IF(ZX) 9995,9995,1000 

9995 READ (5, 11) IA , JA ,IB , JB , SFOC 1 , SFCC2 , V AL1 , VAL2 , J1 , J2 
WRITE (6,11) IA, JA,IE, JB , S FOC 1 , S FOC 2 , VAL 1 , VAL2, J 1 , J2 

1C00 A1 (1, 1)=1.0 
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A1 (1,2) =SF0C1 
A1(2,2) =1.0 
A1 (3, 3) =1.0 
A 1 (4, 4) =1.0 
A 1 (5,5) =1 .0 
A 1 (5, 6) =SFOC1 
A 1 (6,6) =1.0 

DO 1001 L1=1 ,6 
DO 1001 L2=1,6 
TINT 1 (L 1 ,12 ,F) =0 . 0 
DO 1001 13=1,6 

TINT1 (LI ,12 , F) =TINT 1 (L 1 , L 2 , F) ♦ A 1 (L 1 , L3) *T (L 3, L2) 
1001 CONTINUE 

WRITE (6,75) F, (TINT 1 ( 1 , N 1 , F) , N 1 = 1 , 6) 

WRITE (6,75) F, (TINT1 (2,N1,F) ,N1=1,6) 

WRITE (6,75) F, (TINT 1 (3 , N 1 , F) , N 1= 1 , 6) 

WRITE (6,75) F, (TINT 1 (4, N 1 , F) ,N1=1,6) 

WRITE (6, 75) F, (TINT 1 (5 , N 1 , F) , N1 = 1 , 6) 

WRITE (6,75) F, (TINT 1 (6 , N 1 , F) , N 1= 1 , 6) 

A2 (1, 1) =1.0 
A2 (1,2) =SFOC2 
A2(2,2) =1.0 
A2 (3,3) =1.0 
A2(4,4) =1.0 
A2 (5,5) =1.0 
A2 (5, 6) =SFOC2 
A2 (6,6) =1.0 

DO 21002 K 1= 1 ,6 
DO 21002 K2= 1 , 6 
TINT2 (K 1 , K2 , F) =0 . 0 
DO 21002 K 3= 1 , 6 

TINT2 (K1,K2,F)=TINT2(K1,K2,F) ♦A2(K1,K3) *T(K3,K2) 
21002 CONTINUE 

WRITE (6,75) F, (TINT2 ( 1 , N 1 , F) , N 1= 1 , 6) 

WRITE (6,75) F, (TINT 2 (2, N 1, F) ,N1=1,6) 

WRITE (6,75) F, (TINT2 (3 , N 1 ,F) ,N1=1,6) 

WRITE (6,75) F, (TINT 2 (4 , N 1 , F) ,N1=1,6) 

WRITE (6,75) F, (TINT2 (5 , N 1 ,F) , N1= 1 , 6) 

WRITE (6,75) F, (TINT 2 (6 , N 1 , F) ,N1=1,€) 

2000 DIFF1 = VAL1-TINT 1 (IA,»1A, 1) 

DIFF2=VAL2-TINT2 (IB,JB,1) 
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7 FORMAT (/, 3H J= , 12 , 3 X , 4HFHO= , IF 1 C . 4 , 3 X ,4HFHI = , F 10 . 4 , 3X , 6 H ALPH A=, 
1F10.4,3X,5HEETA=,F10.4) 

8 FORMAT (/,3H J=,I2,3X,7H ERIFT=, F10. 4) 

9 FORMAT (/, 3H J= , I 2 , 3X , 1 7HEFFECTI VE LENGTH= , F 1 0. 4 , 5X , 2HK=, F 10 . 4) 

10 FORMAT (10X,6E15.6) 

11 F0RMAT(4I3,4E15.6,2I3) 

12 FORMAT ( 10X,8 H G ( 1 , 1) = , IPS 1 5. 6 , 5X, 8H G ( 1, 2) =, 1PE15. 6) 

13 FORMAT (10X,8H G (2 , 1 ) = , 1 PE 1 5 . 6 , 5 X , 8H G (2 ,2) = , 1PE 15. 6) 

14 FORMAT ( 1H1, IX , 2H J= , 12 , 5X , 5 HX (J) = ,I1) 

15 FORMAT (/,3H J= , 12 , 3X , 1 1 HFOT AT IC N OF, 1F8 . 3, 7HDEGREES) 

16 FORMAT(5F10.5,I1) 

75 FORMAT (1X,2Hf=, 1I2,6E15.6) 

END 


SIEFTC QUAD6 DECK 

SUBROUTINE QUAD (J, EL ,FG) 

DIMENSION A (6,6, 50) ,T (6,6) ,XMAX (6) 

COMMON A ,T , X H AX 
IF(FG) 100,200,300 

100 FG=-FG 

A (1, 1,J) =0.50* (EXP (FG*EI) + EXP (-FG*EL) ) 
A(1,2,J)=(EXF(FG*EL) -EXF (- EL*f G) ) / (2 . 0* FG) 
A (1 , 3 , J) =0 . 0 
A (2, 1, J) =FG*FG*A (1 ,2, J) 

A(2,2,J)=A (1,1, J) 

A (2,3, J) =0.0 
A (3, 3, J) =1.0 
A(4,4,J)=1.0 
A (5,5, J)=COS (FG*EL) 

A (5,6, J) =SIN (FG*EL) /FG 
A(6,5,J)=-FG*SIN (FG*EL) 

A (6,6, J) =A (5,5, J) 

FG=-FG 
GO TO 400 

300 A (1, 1 , J) =CCS (FG*EL) 

A (1,2,0) =SIN (EL*FG) /FG 
A (2, 1 ,J) =-FG*SIN (FG*EL) 

A (2,2,J)=A (1,1, J) 

A (3,3, J) =1.0 
A (4,4, J)=1 .0 

A (5,5, J) =0.50* (EXP (EL*FG) + EXF (-EI*FG) ) 

A (5,6,0) = (EXE (FG*EL) -EXP (-FG*EI) ) / (2.0*FG) 
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A (6,5,J) =FG*FG*A (5,6,J) 
A (6,6, J)=A(5,5,J) 

GO TO 400 
200 A (1,1, J ) = 1.0 
A ( 1 , 2 , J) = EL 
A(1,3,J)=0.0 
A (2, 2, J) =1 . 0 
A (3, 3, J) =1.0 
A(4,4,J) =1.0 
A (5,5, J) =1 .0 
A (5, 6, J) =EL 
A (6 ,6 , J) = 1 . 0 
GO TO 400 
400 RETURN 
END 


SIBFTC EXITU DECK 

SUBROUTINE EXIT (CEL) 

COMMON A ,T, XMAX 

DIMENSION A (6,6,50) ,T (6,6) ,TLAST (6,6) ,AL (6,6) , XMAX (6) ,SIGM (6,6) , 
1 SIG1 (6,6) ,XMF (6) ,TLT (6,6) 

WRITE (6,1) CEL 
DO 200 N= 1 , 25 
Y = N-1 
XL=Y*DEL 


AL(1,1)=1.0 
AL ( 1 , 2) =XL 
AL (2,2) =1.0 
AL ( 3, 3) = 1 . 0 
AL (4,4) = 1.0 
AL (5, 5) =1.0 
AL (5,6) =XL 
AL (6, 6) =1.0 
DO 101 11=1,6 
DO 101 12=1,6 
TLAST(I1,I2)=0.0 
DO 101 13=1,6 

TL AST (II ,12) =T LAST (1 1 , 1 2) ♦ AL ( 1 1 , 13) *T (I 3 , 1 2) 
101 CONTINUE 
WRITE (6,3) 

WRITE (6, 4) XL, (TLAST (1 ,K) ,K=1,6) 

WRITE (6 ,5) (TLAST (2 , K) ,K=1,6) 
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C (M1,M2)=A (M1,M3,JJ)*T (M3, M2) +C (Ml, *2) 

100 CONTINOE 

DO 101 Ml=1,6 
DO 101 M 2= 1 , 6 
T (Ml, M2) =C (Ml, M2) 

101 CONTINUE 
RETURN . 

END 


SIEFTC MAG6 DECK 

SUBROUTINE M AG ( 3 , RO , PHI , A I P , B ET) 

COMMON A,T, XMAX 

DIMENSION A (6,6,50) , T (6 , 6 ) , X M A X (6) 

PHIR=PHI*3. 14 159265/180. 

ALPR=ALP* 3. 14159265/16C. 

BETR=BET*3. 14159265/180. 

A ( 1 , 1 , J) =COS (PHIR-ALPR) /COS (AIFR) 

A(1,2,J)=RC*SIN (FHTR) 

A ( 1, 3, 3) =RO* (1. 0-COS (PHIR) ) 

A (2,1,J)=-SIN (PHIR-ALPR-BETB) / (RO*COS (ALPR) *COS(BETR) ) 

A (2,2,0) =COS (PHIR -BE TR) /COS (BET R) 

A (2,3, J) =SIN (PHIR) ♦ (1 .O-COS (FHIR) ) + T AN (BETR) 

A (3,3, J) = 1 . 0 

A (4, 1 , J) = SIN (PHIR) ♦ (1 . v-CCS (PHIR) ) *T AN (ALPR) 

A (4,2, J)=RC* (1.0 -COS (FHIR) ) 

A (4 , 3 , J) =RO» ( PHIR -SIN (PHIR) ) 

A (4,4, J) =1 .0 

A (5,5,J) =1.C-PHIR*TAN (AIPR) 

A (5,6, J) =FHIF*FO 

A (6, 5, J) = (-TAN (ALPR) -TAN (SETR) ♦FHIR+TAN (ALPR) *T A N (BETR) ) /RO 
A (6,6, J) =1 .0-FHIR*TAN (EETR) 

RETURN 

END 
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Figure 3. - Evolution of one-dimensionalphase-space ellipse as beam envelope passes through a 
waist Lower diagram indicates physical size of beam envelope over corresponding interval. 



Figure 4 - Schematic of PDP-15 program flow. 
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Figure 6. - Typical PDP-15 display of output from BEAMR. System features two sets of quadrupole doublets 
with a dipole magnet between them. 
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